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Abstract. 

We formulate and solve a two-state model for the elasticity of nicked, double-stranded DNA 
that borrows features from both the Worm Like Chain and the Bragg-Zimm model. Our model 
is computationally simple, and gives an excellent fit to recent experimental data through the 
entire overstretching transition. The fit gives the first value for the bending stiffness of the 
overstretched state as about 10 nm ■ feT, a value quite different from either B-form or single- 
stranded DNA. 



Introduction and Summary. - When double-stranded DNA is subjected to longitudinal 
forces greater than about 65 pN it undergoes a radical conformational change, marked by a 
sudden, almost twofold increase in contour length G],[|. The structural characterization of the 
resulting overstretched or ("S-state") DNA is complicated by the fact that techniques such 
as X-ray crystallography are not applicable to single molecules. In this Letter, we instead 
characterize overstretched DNA by examining its elastic constants, and to this end formulate 
and solve a model that synthesizes features of both the Worm Like Chain (WLC) and the 
Bragg-Zimm model of the helix-coil transition in peptides. Thus we model DNA as consisting 
of two different, coexisting conformations, each with its own elastic constants. We solve this 
model and show that it gives a good fit to recent data on the overstretching transition in 
nicked, double-stranded DNA. From these fits, we conclude that the bend stiffness of S-DNA 
is intermediate between the known values for single stranded and double stranded DNA. Our 
result supports the work of Leger et al. ||, || , who argued that S-DNA has a definite helical 
pitch and hence is a new duplex conformation of DNA. 

Our model and solution method differ from those offered by Marko who assumes the 
bend stiffnesses of the two conformational states to be identical; our analysis will show that 
on the contrary the stiffnesses are markedly different. The analysis of Viovy and Cizeau [|| is 
essentially a mean-field approximation to the model we study here; in addition, the authors 
did not quote any value for the S-DNA bend stiffness, presumably because the experimental 
data available at that time did not permit such a determination. 
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The model studied here is a continuum limit of a more general class of discrete persistent- 
chain models. Such models give better fits to the stretching curves of smg/e-stranded DNA 
than either the continuum WLC or the freely jointed chain models. Details will appear 
elsewhere 0. 

Our model and method are also of some general interest beyond DNA. For example, both 
can be adapted to the study of the stretching of polypeptides with a helix-coil transition. 

Model. - We begin by formulating a discretized form of the WLC, which we call the the 
"Discrete Persistent Chain." Later we will introduce an Ising-like variable on each chain link 
describing a cooperative transition from B- to S-form. 

The DPC models the polymer as a chain of N segments of length b, whose conformation 
is fully described by the collection of orientation vectors {ti} for each segment. Thus the 
relaxed total contour length is L tot = Nb. Bend resistance is taken into account by including 
an energy penalty at each link proportional to the square of the angle S^j+i = arccos(ij • ij+i) 
between two adjacent links. The energy functional describing this model is thus given by 

v fb i 2 + V A (q ) 2 m 

i—1 i=l 



s £[{t,}]/fc B T where §2 



The partition function for this energy functional is Z 
is the two-dimensional unit sphere. 

To compute Z we use the transfer matrix technique J8j, interpreting each integral above 
as a generalized matrix product among matrices with continuous indices: Z — v ■ T^ -1 • w . 
In this formula v and w are vectors indexed by t, or in other words functions v(t), w(t ). The 
matrix product T • v is a new vector, defined by the convolution: 



(T.tf)(t 4 )= / d^T^tX^). (2) 

Js 2 

Here the matrix elements of T are given by J(ti,ij) = e~ £i ^ ti,t i^ knT ; we will not need the 
explicit forms of v and w below. 

The force-extension relation can be obtained from Z by differentiating with respect to 
the force. It is here that the transfer matrix formulation can be used to greatly simplify the 
calculation of the force-extension relation, since all that is needed to compute the logarithmic 
derivative of Z in the limit of long chains is the largest eigenvalue of T, which we will call 

Amax- 

Z v largeJV fk B T\ d , N ( k B T\ d 



Analogously to ref. Q, it is straightforward to add an intrinsic stretch modulus to the 
calculation outlined above, obtaining an "Extensible DPC" model. 

To study overstretching, we now extend the extensible DPC by giving each link a discrete 
variable a, which takes the values ±1. We will take a = +1 to mean the segment is in the 
B-state and a = — 1 for the S-state. The factor by which a segment elongates when going from 
B to S will be called C, i-e. 6 (S) = £b (with C > 1). We assign a bend stiffness parameter A to 
B-DNA, and a different = f3(A to S-DNA; (3 is a dimensionless parameter with /3£ < 1. 
Similarly we assign a bend stiffness r\A to a hinge joining a B and an S segment. 
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The full energy functional for the Ising-DPC model is thus: 
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The first line is the pure-Ising part, with 2aok B T the intrinsic free energy cost of converting 
a single segment from B to S and 2 r yk B T the energy cost of creating a B— >S interface. Note 
that we ignore a contribution to the energy functional from the first and last segments. In 
the long-chain limit this does not affect the outcome of our calculation. 
The partition function for the energy functional (M) is given by 



N-l 

n e 
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where now £[{tj, cr^}] = X^ili 1 £j(*i> ViiU+i, Cj+i)- We again calculate Z with the aid of the 
transfer matrix technique, writing eq. ^| as Z = v ■ T N ~ 1 • w , with T now the transfer matrix for 
our Ising-DPC model, which carries an additional 2-by-2 structure due to the Ising variables. 
The dot products are thus defined as 



(T-fl)«r 4 (&)= 52 / d%- T ffiCTj (M> CTj (<;■)* 



(6) 



<7j=±l 

The individual matrix elements T aiCrj are given explicitly by 

1 



T_i,-i(ii, U+ij 



exp 



(3A 

I 



-a 



and related expressions for Tx,i, Ti._i, and T_x,i, where / = j^pf- 

We approximate the largest eigenvalue of the transfer matrix T using a variational approach 
flTof . We choose a three- parameter family of trial eigenfunctions with azimuthal symmetry, 
peaked in the direction of the force z: 



( Ul V 

\^sinh(2wi) J 

( »-i v 

\^sinh(2w_i) J 



(7) 



These trial functions were chosen such that their squared norm is independent of all pa- 
rameters: ||?V,w_i,y|| 2 = 2n. Eq. ^ shows that the o/s give the degree of alignment of the 
monomers (how forward-peaked their probability distribution is), whereas ip describes the 
relative probability of a monomer to be in the two states. The variational estimate for the 
maximal eigenvalue is thus 



A max= max y(ul,Ul_x,lp) 

Wl ,CJ_ 1 ,</? 



max 
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(8) 
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The maximization over ip can be done analytically: defining the 2x2 matrix T(a>i,a;_i) 

(9) 



by 

gives that 



Vuttu-uv ' T ' ^wi,w_i,v> = (cos^, sinyj) • T(wi,w_i) 



cos^ 
simp 



A max = max 



12 ' 



(10) 



where y(wi,w_i) is the maximal eigenvalue of the 2x2 matrix T(wi,w_i). The following 
section will calculate this eigenvalue in a continuum approximation to T(u>i, illustrating 
the procedure by considering in some detail the matrix element Ti t x(u>i,u>—i). The other 
matrix elements can be obtained analogously. Writing out the integrals explicitly, we have 



sinh(2u;i) 



s 2 



d 2 t 



i+l 



(11) 



where we have introduced a 
f) 2 + 2fif 



u 1 



Condensing notation even further we define /i 2 



a 2 + (4) 2 + 2a4fi • z, which allows us to write 



ti,i(wi) = (2 



sinh(2 Wl ) /Aju aA 



(12) 



Continuum Limit. We could now proceed to evaluate the force-extension relation 
of the Ising-DPC model, by evaluating eq. [l^ numerically and using eq. |^. To simplify the 
calculations, however, we first pass to a continuum limit. To justify this step, we note that the 
continuum (WLC) approximation gives a good account of single-stranded DNA stretching out 
to forces beyond those probed in overstretching experiments (about 90 pN) jflj . As mentioned 
earlier, the continuum approximation is also quite good for double-stranded DNA, because 
the latter's persistence length is so much longer than its monomer size. 

In the continuum limit b is sent to zero holding L to t fixed; hence N — > oo. The bookkeeping 
is more manageable after a shift in /i: x = /i — (A/b). Eq. Ilj then reduces to 
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The last integral can be worked out exactly, and expanding the result to second order in b we 
end up with 
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In similar fashion, we can obtain the following expressions for the other matrix elements. 
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To obtain a nontrivial continuum limit we must now specify how the parameters A, ao, 
and 7 depend on b as b — > 0. The choices 

a Q = --lnfl + ba, 7 = --ln(g6) (16) 

give a well-defined limit, where we hold A, a, (3 and g fixed as b — ► 0. With these choices, the 
matrix m rr^T(u)i.u>—i) takes the form 



with 



g = ffV? / o^-i \ 5 / 2sinh(ui +cj_i) , 



sinh(2wi) sinh(2o;_i)/ \ lu± + ui-\ 

Note that the prefactor -£^=j in eq. [Tl] does not contribute to the force-extension result eq. ||, 
because it does not depend on the force. In terms of the individual matrix entries, the quantity 
to be maximized now reads (see eq. ^|): 

in jr(«i, w_o = h - {v + n + VCP - ny + 4Q 2 ) . (19) 

Writing il = In A* nax = b^ 1 x maxlny(wi, a;_i), the force-extension in the continuum limit 
is finally given by 

(-T-) = k B T§ . (20) 

We evaluate Q by numerically maximizing eq. [l9| 

So far, we have not included stretch moduli for the B- and S-DNA. This is easily imple- 
mented to first order in f / E by replacing / with /(l + 2E { s ,b) ) in the matrix elements for the 
two states respectively (eq. 7). This procedure yields theoretical force-extension curves like 
the one plotted in Fig. ([y). 

In summary, our model contains the following seven parameters. 2ak B T is the free energy 
per unit length required to flip B-DNA into the S-state, and is measured in [J/nm]. Q measures 
the cooperativity of the transition and has units [1 /nm] . A is the bend stiffness parameter of 
B-DNA, with units [nm]. The dimensionless parameter f3 is the ratio of the B- and S-DNA 
bend stiffnesses. E^ and E^ are the stretch stiffnesses of B and S-DNA, and are measured 
in pN. Finally, £ is the dimensionless elongation factor associated with the B— >S transition. 

Discussion of fit. - Fig. (|j) shows a fit to some recent experimental data (similar data 
appear in E^|). Our model reproduces the experimental data rather well, but with so many 
fit parameters one may ask whether it actually makes any falsifiable predictions. To answer 
this question we note that the data below the transition suffice to fix A and E^ as usual, 
roughly speaking from the curvature and slope of the curve below the transition. Similarly, 
the data above the transition fix A^ — C,f3A and E^ s \ The vertical jump in the curve at 
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the transition fixes £. The horizontal location of the jump fixes a, and the steepness of the 
jump fixes the cooperativity Q.(^j Thus all of the model's parameters are fixed by specific 



features of the data. Two additional, independent features of the data now remain, namely 
the rounding of the curve at the start and end of the transition. Our model predicts these 
features fairly succesfully. 

The fit recovers the known values for the effective persistence length of B-DNA of around 
50 nm and its stretch modulus of about fOOOpN. Our first result is that the bend stiffness of 
S-DNA from our fit as A^ = = 12.32 nm. Similar results were obtained using the older 
data of Cluzel et al. |l],|7]]. If S-DNA consisted of two unbound, single strands, we might have 
expected to be twice as large as the value A ss w 0.75 nm appropriate to single-stranded 
DNA (as obtained from stretching experiments on ssDNA, restricted to forces above those 
required to pull out secondary structure [^} |ll],|7) ) . On the contrary, we find that the bend 



stiffness of S-DNA is intermediate between that of B-DNA and that of two single strands. 

Finally, our fit gives the stretch modulus of S-DNA is substantially higher than that of 
B-DNA. This conclusion is consistent with the idea that the contour length of S-DNA is 
determined by its covalently bonded sugar-phosphate backbones, which are much straighter 
than in B-DNA; the contour length of B-DNA is instead determined by weaker, base-stacking 
interactions. 

Relation to prior work. - Several authors have studied the entropic elasticity of two-state 
chains. As soon as the overstretching transition was discovered, Cluzel proposed a pure Ising 
model by analogy to the helix-coil transition |lj] . Others then introduced entropic elasticity, 
but required that both states have the same bending stiffness as B-DNA ||[l5) or took one of 



the two states to be infinitely stiff 16 1, or to be a FJC [17]]. Also several earlier works made 
a mean-field approximation instead of diagonalizing the full transfer matrix. We believe our 
Ising-DPC model to be the first consistent formulation incorporating the coexistence of two 
different states with arbitrary elastic constants. Our approach also is calculationally more 
straightforward than some, and minimal in the sense that no unknown potential function 
needs to be chosen. 

* * * 

We thank T. Burkhardt, D. Chatenay, A. Grosberg, R. Kamien, J. Marko and M. Rief for 
valuable discussions, and C. Bustamante, D. Chatenay, J.-F. Leger, J. Marko, M. Rief, and 
S. Smith for sending us experimental data. CS acknowledges support from NIH grant R01 
HL67286 and from NSF grant DMR00-79909. PN acknowledges partial support from NSF 
grant DMR98-07156. 



REFERENCES 

[1] Cluzel P., Lebrun A., Heller C, Lavery R., Viovy J.-L., Chatenay D. and Caron F., 

Science, 271 (1996) 792. 
[2] Smith S. B., Cui Y. and Bustamante C, Science, 271 (1996) 795. 

[3] Leger J.-F., Romano C, Sarkar A., Robert J., Bourdieu L., Chatenay D. and Marko 
J. F., Phys. Rev. Lett, 83 (1999) 1066. 



( 1 ) The fit value of a should be regarded as an average of the two different costs to convert AT or GC pairs. 
The fit value of Q has no direct microscopic significance, as the apparent cooperativity of the transition will 
be reduced by the sequence disorder. 

( 2 ) Hagerman's result that the persistence length of a single strand of poly(dT) DNA is between 2 nm and 3 
nm pjj does not come from a stretching experiment and should not be compared directly. 



C. Storm and P. C. Nelson: The bend stiffness of S-DNA 



7 



[4] Leger J.-F., L'ADN : une flexibthte structurale adaptee aux interactions avec les autres macro- 
molecules de son environnement (Ph.D. thesis, Universite Louis Pasteur, Strasbourg) 1999. 
[5] Marko J. F., Phys. Rev. E, 57 (1998) 2134. 
[6] Cizeau P. and VlOVY J.-L., Biopolymers, 42 (1997) 383. 
[7] Storm C. and Nelson P. C.,. In preparation, (2002) . 
[8] Kramers H. A. and Wannier G. H., Phys. Rev., 60 (1941) 252. 
[9] Odijk T., Macromolecules, 28 (1995) 7016. 

[10] Marko J. F. and Siggia E. D., Macromolecules, 28 (1995) 8759. 

[11] Clausen-Schaumann H., Rief M., Tolksdorf C. and Gaub H. E., Biophys. J., 2000 (78) 
1997. 

[12] Bustamante C, Smith S. B., Liphardt J. and Smith D., Curr. Op. Str. Biol, 10 (2000) 
279. 

[13] Mills J. B., Vacano E. and Hagerman P. J., J. Mol. Biol, 285 (1999) 245. 

[14] Cluzel, P., L'ADN, une molecule extensible (Ph.D. thesis, Universite Paris VI) 1996. 

[15] Ahsan A., Rudnick J. and Bruinsma R., Biophys. J., 74 (1998) 132. 

[16] Tamashiro M. N. and PlNCUS P., Phys. Rev. E, 63 (2001) 021909. 

[17] Rouzina I. and Bloomfield V. A., Biophys. J., 80 (2001) 882. 




0.5 



I i i i i I 

20 40 60 80 100 
f [pN] 

Fig. 1 - Least-squares fit of the Ising-DPC model to an overstretching dataset (48.5 kbp A DNA 
construct; buffer 500 mM NaCl, 20mM Tris, pH 8). Data kindly supplied by C. Bustamante and 
S. Smith. Fit parameters: A = 43.75 nm, a = 5.45 nm -1 , j3 = 0.16, Q = 0.13 nm" 1 , £ = 1.76, 
£ (B) = 1.2 ■ 10 3 pN and E (s) = 1.0 ■ 10 4 pN. x = 9.22 at N = 825; points with 1.11 < (f ) < 1.55 
were excluded from the fit. 



